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To  predict  the  warpage  behaviour  of  an  anode-supported  half-cell  of  a  solid  oxide  fuel  cell  (SOFC),  a 
viscoelastic  analysis  has  been  made  for  a  system  with  Gd-doped  ceria  (GDC)  and  NiO-yttria  stabilized 
zirconia  (NiO-YSZ)  bi-layers.  The  viscoelastic  properties  of  each  component  at  1573  K  are  measured  by 
discontinuous  sinter  forging  (for  GDC)  and  continuous  sinter  forging  (for  NiO-YSZ)  techniques.  Using 
the  measured  uniaxial  viscosities  and  the  viscous  Poisson  ratio,  the  warpage  during  the  co-firing  of  a 
GDC/NiO-YSZ  bi-layer  is  calculated  for  the  continuum  mechanical  models  of  Cai  et  al.  and  Ranters  et  al. 
The  warpage  predicted  by  the  two  models,  in  particular  that  of  Cai  et  al.,  is  shown  to  be  in  good  quantitative 
agreement  with  the  experimental  observation.  This  result  demonstrates  that  the  warpage  of  bi-layer 
samples  during  co-firing  can  be  suitably  predicted  when  the  viscoelastic  properties  of  each  component 
are  properly  measured  and  estimated.  The  present  investigation  also  provides  a  useful  methodology  for 
optimizing  the  design  of  planar  asymmetric  half-cells  for  the  SOFC  industry. 

©  2008  Elsevier  B.V.  All  rights  reserved. 


1.  Introduction 

A  planar  geometry  of  stacked  components  is  widely  chosen  for 
use  in  solid  oxide  fuel  cells  (SOFCs)  due  to  its  high  power  density  and 
low  fabrication  cost  [1,2].  One  of  the  critical  issues  in  fabricating  a 
planar  anode-supported  SOFC  is  the  co-firing  of  the  solid  electrolyte 
and  anode.  During  the  co-firing  of  asymmetric  bi-layers,  sintering 
defects  such  as  cracks,  debonding  and  warpage  often  occur  due 
to  the  different  densification  rates  that  exist  among  the  compo¬ 
nents.  In  order  to  solve  this  problem,  the  sintering  behaviour  of  each 
component  has  to  be  characterized  [3-11].  Nevertheless,  few  con¬ 
tinuum  mechanical  approaches  have  been  attempted  in  an  effort 
to  understand  the  co-firing  processes  of  SOFC  components. 

Another  issue  in  the  fabrication  of  planar  anode-supported 
SOFCs  is  related  to  their  operating  temperature.  As  the  operation 
temperature  of  the  well-known  zirconia-based  electrolyte  is  too 
high  (~1273  K)  for  continuous  use  for  a  long  period  of  time,  other 
electrolyte  materials  have  been  investigated  [1].  Gd-  or  Sm-doped 
ceria  is  each  considered  to  be  a  good  candidate  for  low-temperature 
operating  SOFCs  on  account  of  the  high  electrical  conductivity  of 
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~0.08Scm-1  at  1073  K  [12,13].  As  an  anode  material,  a  NiO-yttria 
stabilized  zirconia  (denoted  as  NiO-YSZ)  composite  is  commonly 
used  because  of  its  good  mechanical  strength  and  chemical  stability 
at  a  low  partial  pressure  of  oxygen  [1,14]. 

To  determine  the  viscoelastic  behaviour  of  a  porous  multi-layer 
body  throughout  a  densification  process,  the  viscoelastic  proper¬ 
ties  of  individual  layers  initially  have  to  be  measured.  Scherer  [3] 
suggested  that  viscoelastic  properties  such  as  the  uniaxial  viscos¬ 
ity  and  viscous  Poisson  ratio  could  be  determined  from  a  linear 
viscous  constitutive  relationship  of  the  mechanical  response  of  a 
porous  body.  The  uniaxial  viscosity  and  viscous  Poisson  ratio  as  a 
function  of  the  porosity  are  obtained  by  means  of  continuous  [4], 
cyclic  [5]  and  discontinuous  sinter  forging  techniques  [6-8]. 

The  purpose  of  the  present  investigation  is  to  predict  the 
warpage  behaviour  of  electrolyte] anode  bi-layer  samples  using  the 
measured  viscoelastic  properties  of  the  components  in  a  continuum 
mechanical  approach  and  to  test  the  applicability  of  the  predic¬ 
tion  by  observing  the  warpage  behaviour  of  bi-layer  samples.  The 
viscoelastic  properties  of  nano-sized  Gd-doped  ceria  (GDC)  and 
NiO-YSZ  are  measured  by  discontinuous  and  continuous  sinter 
forging  techniques,  respectively,  and  the  warpage  behaviour  of  the 
GDC/NiO-YSZ  bi-layers  during  co-firing  is  evaluated.  The  observed 
warpage  of  the  bi-layers  at  1573  K  is  compared  with  calculations 
using  models  proposed  by  Cai  et  al.  [9]  and  Ranters  et  al.  [15]. 
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2.  Theoretical  background 


3.  Experimental  procedure 


In  cylindrical  coordinates,  the  axial  and  radial  strain  rates  of  a 
sintering  body,  ez  and  £r,  are  expressed,  respectively,  as  [3]: 

£z  =  Sf  +  Ep1  [<TZ  —  Vp(crr  +  CTq)]  (1) 

and 

£r  =  +  Ep^[crr  —  Vp(crQ  +  crz)]  (2) 

Assuming  that  the  hoop  and  radial  stresses,  cre  and  crr,  are  negligible 
under  uniaxial  loading,  the  uniaxial  viscosity,  Ep  is  defined  by: 


where  oz  is  the  applied  stress  and  £f  is  the  free  strain  rate.  Under 
the  same  assumption,  the  viscous  Poisson  ratio,  vp,  is  expressed  as 
[3]: 


vp  - 


gr-gf 
gZ  -gf 


(4) 


A  more  reliable  method  for  determining  vp  has  been  derived  from 
the  volumetric  densification  rate,  £v,  under  uniaxial  loading  [16]: 


£y  —  £z  +  2£f  —  3Sf  +  -=r~  (1  —  2vp) 
Ep 


(5) 


Given  Ep,  the  viscous  Poisson  ratio  can  be  determined. 

Warping  of  an  asymmetric  bi-layer  is  due  to  the  difference  in  the 
shrinkage  rate  between  the  two  layers,  which  results  in  mismatch 
stresses  [9].  Different  continuum  mechanical  models  can  be  used 
to  predict  the  camber,  all  based  on  Kirchhoff  assumptions  (nor¬ 
mals  remain  straight,  unstretched  and  normal)  for  a  thin  plate  and 
neglecting  edge  effects.  The  radial  strain  of  the  laminate  £'r  can  be 
defined  as  [15]: 


£'r(z)  =  £q  -  zl<! 


(6) 


where  k  is  the  laminate  camber  (inverse  of  the  curvature  radius) 
and  £q  is  the  strain  for  z  =  0  (by  convention  laminate  mid-plane). 
Combining  equations  describing  strain  rates,  Eqs.  (1)  and  (2),  the 
stress  along  the  laminate  thickness  is  expressed  as: 

<xr(z)  =  ,  (*0  -  2k!  -  zK  -  £f(z))  (7) 

1  —  Vp(Z) 


By  applying  the  force  and  torque  equilibrium  on  the  entire  system 
(discretized  in  500  sub-layers),  the  two  unknowns,  \<!  and  £q,  are 
deduced.  The  radial  and  axial  strain  rates  are  computed  from  Eqs. 
(1)  and  (2)  and  from  the  density  distribution  along  the  thickness 
axis.  Sintering  parameters  (Ep,  vp  and  £f)  are  then  corrected  with 
the  new  calculated  density  for  the  next  calculation  step. 

By  neglecting  the  term  zl<!  in  Eq.  (7)  and  the  built-in  stress 
that  affects  densification,  the  model  of  Cai  et  al.  [9]  allows  a  more 
straightforward  calculation.  This  is  given  by: 


k  = 


d((ti  +t2)/r) 
d  t 


6 (m  -  1  )2mn 

m4n2  +  2mn(2m2  +  3m  +  2)  +  1 


A£ 


(8) 


3.1  Sample  preparation  for  viscous  property  measurement 

3.1.1.  Gd-doped  ceria 

Nanocrystalline  Gd-doped  ceria  (Ce0.9Gd0.iO2_5,  denoted  as 
GDC,  Rhodia,  Germany)  powder  was  used.  The  average  size  of  the 
powder,  as  measured  by  the  BET  method  (Autosorb  3B,  Quan- 
tachrome),  was  32  nm.  As  this  Gd-doped  ceria  powder  is  nano-sized 
with  a  high  driving  force  for  densification,  the  densification  was 
mostly  completed  upon  reaching  the  co-firing  temperature  of 
1573  K.  Thus,  the  viscoelastic  properties  were  measured  at  lower 
temperatures  and  they  were  estimated  at  higher  temperatures 
through  the  use  of  appropriate  functions  of  density,  grain  size  and 
activation  energy.  Using  a  12  mm  diameter  stainless-steel  die,  cylin¬ 
drical  pellets  were  uniaxially  dry-pressed  at  100  MPa.  After  cold 
isostatic  pressing  at  350  MPa  for  120  s,  the  green  body  density  was 
60.7  ±  0.3%  of  the  theoretical  value. 

Cylindrical  pellets  of  a  height  of  16  mm  were  employed  in  sinter 
forging  experiments  using  high-resolution  axial  and  radial  lasers  to 
measure  the  height  and  diameter  of  the  samples  [17].  Some  pellets 
were  freely  sintered,  whereas  others  were  sintered  up  to  a  desired 
density  at  1323, 1373  and  1423  K.  They  were  then  temporarily  sub¬ 
jected  to  a  constant  load  of  either  50  or  100  N,  corresponding  to 
0.6  and  1.2  MPa,  respectively.  This  mechanical  load  was  applied 
within  a  5%  increase  of  relative  density.  The  heating  and  cooling 
rate  was  set  at  30  K  min-1 .  The  density  of  each  sintered  sample  was 
determined  by  the  Archimedes  method.  The  microstructure  was 
characterized  on  fractured  surfaces.  The  grain  size  was  measured 
by  the  lineal  intercept  method  with  a  correction  coefficient  of  1.56 
using  LINCE  software  (ver.  2.31,  Ceramics  Group,  TU  Darmstadt).  As 
the  mean  grain  size  obtained  from  a  fractured  surface  is  smaller 
than  that  measured  on  a  plane  section  by  ~7%  [15],  compensation 
was  made  when  presenting  the  data.  More  than  150  grains  were 
examined  for  each  sample. 

3.1.2.  NiO-YSZ  composite 

NiO  (Sumitomo,  Japan)  and  YSZ  (TZ-8Y,  Tosoh,  Japan)  were  used 
to  prepare  a  composition  of  40vol.%  of  NiO  and  60vol.%  of  YSZ. 
The  mixed  powders  were  ball-milled  in  ethanol  for  24  h  using 
5  mm  diameter  zirconia  balls.  Green  bodies  12  mm  in  diameter 
and  20  mm  in  height  were  compacted  using  a  stainless-steel  die 
at  100  MPa  and  were  then  cold  isostatically  pressed  at  175  MPa  for 
120  s.  The  green  density  of  the  NiO-YSZ  composite  was  52.3%  of  the 
theoretical  value  determined  by  the  rule  of  mixture.  The  shaped 
pellets  were  freely  sintered  at  1573  K,  and  some  of  them  were  then 
continuously  subjected  to  a  constant  load  of  100  N,  corresponding 
to  ~1.2  MPa.  As  the  continuous  and  discontinuous  sinter  forging 
techniques  gave  similar  results  for  alumina  samples  with  a  relative 
density  lower  than  80%  [18],  a  continuous  sinter  forging  technique 
was  adopted.  This  technique  requires  fewer  specimens  compared 
with  a  discontinuous  sinter  forging  technique.  The  heating  and 
cooling  rate  was  set  to  30  K  min-1.  The  density  of  each  sintered 
specimen  was  determined  by  the  Archimedes  method. 

3.2.  Observation  ofwarpage  behaviour  of  a  GDC/NiO-YSZ  bi-layer 


where  k  is  the  rate  of  the  curvature  change  normalized  to  the  total 
thickness,  Ac  is  the  mismatch  between  the  free  strain  rates  of  two 
layers,  m  is  the  layer  thickness  ratio,  t\  and  t2  are  the  thickness 
of  the  two  layers,  r  is  the  curvature  radius  of  the  bi-layer,  n  is  the 
viscosity  ratio  between  the  layers  and  is  expressed  as  [9]: 


Epi  1  —  vp2 

1  —  vpi  £p2 


(9) 


Bi-layers  of  GDC  and  NiO-YSZ  composite  were  prepared  by 
a  screen-printing  method.  The  NiO-YSZ  composite  powder  was 
dry-pressed  at  100  MPa  for  5  min  in  a  3  mm  x  30  mm  rectangular 
stainless-steel  mold  and  pre-sintered  at  1273 1<  for  30  min  to  obtain 
a  relative  density  of  52.3  ±  0.2%.  The  ethanol-based  slurry  of  GDC 
was  then  screen-printed  on  to  the  prepared  substrate  composites. 

The  bi-layer  was  sintered  at  1573  K  in  a  tubular  furnace  and  its 
shape  was  photographed  through  a  transparent  window  with  a  CCD 
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Fig.  1.  Densification  curves  of  GDC  at  three  different  temperatures. 

camera  every  15  s.  The  heating  rate  was  30  K  min-1 .  The  radii  of  a  bi¬ 
layer  were  analyzed  by  an  image  analysis  program  (ImageJ,  ver.  1.37, 
National  Institute  of  Health).  The  thickness  changes  of  the  two  lay¬ 
ers  with  sintering  time  were  measured  by  contact  dilatometry  (DIL 
402C,  Netzsch,  Germany)  using  1.0  mm  x  1.0  mm  square  samples 
and  identical  heating  schedules.  (At  the  beginning  of  isothermal 
sintering,  the  thickness  of  the  anode  and  the  electrolyte  were 
1.8  mm  and  24  fxm,  respectively.)  With  these  small-size  samples,  it 
was  possible  to  minimize  the  deflection  due  to  the  difference  in  the 
axial  shrinkage  between  the  two  layers.  The  measured  curvature 
was  normalized  to  the  measured  thickness. 

4.  Results  and  discussion 

4. 1.  Viscous  properties  of  Gd-doped  ceria 

Densification  curves  of  the  GDC  pellets  are  shown  in  Fig.  1.  The 
sintering  temperature  could  be  reduced  from  the  conventional  sin¬ 
tering  temperature  of  1773-1973  to  1373  K  using  nano-sized  GDC 
powder,  as  in  the  case  of  Li  et  al.  [19].  The  sintered  density  can  be 
determined  at  any  time  by  the  following  relationship: 


Relative  Density  (%) 

Fig.  2.  Overall  uniaxial  viscosities  with  relative  density  of  GDC  at  different  temper¬ 
atures. 


Relative  Density  (%) 

Fig.  3.  Sintering  trajectories  of  GDC  at  three  different  temperatures. 

Here,  p0  is  the  initial  density  and  sz  and  er  are  the  true  strains.  The 
relative  densities  of  the  samples  sintered  at  1323, 1373  and  1423  K 
for  2  h  were  94.5,  97.3  and  98.4%,  respectively. 

The  uniaxial  viscosity  was  determined  from  the  linear  relation¬ 
ship  between  the  applied  load  and  the  axial  strain  rate.  The  overall 
uniaxial  viscosities  as  a  function  of  the  relative  density  at  three 
different  temperatures  are  shown  in  Fig.  2.  The  uniaxial  viscosity 
increases  as  the  sintering  proceeds  and  drastically  increases  when 
the  specimen  approaches  its  final  density  at  each  temperature.  The 
drastic  increase  in  the  uniaxial  viscosity  is  related  to  grain  growth 
at  the  late  stage  of  sintering  [8]. 

The  uniaxial  viscosity  can  be  expressed  as  independent  func¬ 
tions  of  the  grain  size  and  relative  density  along  with  a  thermally 
activated  term  [8,20,21  ]: 

Ep=aEp(p)Ep(d)e^kT)  (11) 

Here,  a  is  a  scaling  factor,  and  Ep{p)  and  Ep(d)  are  functions  of  the 
relative  density  and  grain  size,  respectively.  For  the  same  Gd-doped 
ceria,  the  function  of  density,  Ep{p),  was  previously  found  to  follow 
Beere’s  model  [8].  The  function  of  the  grain  size,  Ep(d),  is  known 
to  be  proportional  to  the  cube  of  the  grain  size  normalized  to  the 
initial  grain  size  for  grain  boundary  diffusion  [8,20,21],  which  is 
the  case  for  this  GDC.  Fig.  3  plots  the  sintering  trajectories  of  sam¬ 
ples  sintered  freely  at  different  temperatures.  During  the  sintering 
process,  the  grain  size  increased  from  150  to  400  nm. 


Fig.  4.  Arrhenius  plot  of  uniaxial  viscosity  normalized  by  grain  size  as  function  of 
inverse  temperature. 
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Fig.  5.  Representative  plot  of  volumetric  densification  rate  as  function  of  applied 
stress  to  determine  viscous  Poisson  ratio  at  1323  K. 


The  uniaxial  viscosity  normalized  to  the  grain  size  was  obtained 
using  Eq.  (11).  This  is  given  in  Fig.  4  as  a  function  of  the  inverse 
temperature.  Its  apparent  activation  energy  is  estimated  to  be 
3.57  ±  0.05  eV  from  the  slope  of  the  Arrhenius  plot. 

The  viscous  Poisson  ratio  has  also  been  determined  by  a  discon¬ 
tinuous  sinter  forging  technique.  Fig.  5  shows  a  linear  dependence 
between  the  applied  stress  and  (ev  -  £f)Fp  at  80,  85  and  90%  of 
the  theoretical  density  at  1323  K.  The  slopes  in  the  plots  represent 
1  -  2vp.  The  measured  variation  of  the  viscous  Poisson  ratio  with 
the  relative  density  at  1323 1<  is  presented  in  Fig.  6.  Each  set  of  data 
points  was  obtained  from  a  sample  with  a  loading  applied  within 
a  5%  increase  of  the  relative  density.  For  the  range  of  80-84%  rela¬ 
tive  density,  the  measured  values  are  in  good  agreement  with  the 
isotropic  Venkatachari  model,  as  has  been  observed  for  alumina 
compacts  [22].  As  the  relative  density  increases  above  this  range, 
the  deviation  from  the  model  becomes  larger.  This  may  be  due  to 
the  increased  microstructural  anisotropy  induced  by  the  uniaxial 
mechanical  loading.  A  significant  increase  of  the  anisotropic  viscous 
Poisson  ratio  has  also  been  estimated  by  means  of  discrete  ele¬ 
ment  modelling  for  uniaxially  loaded  alumina  compacts  (direction 
3  being  the  direction  of  the  applied  stress)  [23]. 


Fig.  7.  Densification  curves  (a)  and  free  strain  rates  (b)  of  GDC  and  NiO-YSZ  com¬ 
posite  at  1573  K. 


4.2.  Viscous  properties  of  NiO-YSZ 

Fig.  7(a)  plots  the  free  densification  curves  of  NiO-YSZ  (and  GDC 
for  comparison)  at  1573  K  as  a  function  of  the  isothermal  sintering 
time.  Fig.  7(b)  shows  the  free  strain  rate  as  a  function  of  the  density 
for  40NiO-60YSZ  together  with  that  for  GDC  at  1573  K.  The  free 
strain  rate  of  NiO-YSZ  is  higher  than  that  of  GDC,  as  the  relative 
density  of  the  composite  sintered  at  1573  K  for  30  min  was  only  82% 
of  the  theoretical  value,  whereas  GDC  becomes  nearly  fully  densi- 
fied  after  200  s.  The  measured  uniaxial  viscosity  of  NiO-YSZ  as  a 
function  of  density  is  given  in  Fig.  8.  The  uniaxial  viscosity  of  GDC 
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Fig.  8.  Uniaxial  viscosities  of  NiO-YSZ  (measured)  and  GDC  (estimated)  as  function 
of  density  at  1573  K. 
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Fig.  9.  Images  of  GDC/NiO-YSZ  bi-layer  during  co-firing. 


calculated  using  Eq.  (11 )  is  also  plotted.  The  uniaxial  viscosity  value 
is  nearly  80GPas-1  at  82%  of  the  theoretical  density.  Given  that 
the  viscous  Poisson  ratio  depends  mainly  on  the  relative  density 
irrespective  of  the  temperature  and  grain  size  [15,21  ],  the  Venkat- 
achari  model  was  also  used  to  predict  the  sintering  behaviour  of 
the  NiO-YSZ  composite. 

4.3.  Evaluation  ofwarpage  behaviour  of  GDC/NiO-YSZ  bi-layer 

Images  of  warped  bi-layers  with  a  GDC  layer  on  the  top  dur¬ 
ing  co-firing  are  presented  in  Fig.  9.  As  the  GDC  layer  sinters  faster 
than  the  NiO-YSZ  layer,  the  sample  cambers  upwards  first.  When 
the  NiO-YSZ  anode  begins  to  sinter  appreciably,  the  curvature 
decreases  and  eventually  the  sample  flattens.  Fig.  10  plots  the  vari¬ 
ation  of  the  curvature  radius  of  a  bi-layer  sample  while  the  sample 
is  heated  to  1573  K  and  upon  isothermal  holding  at  1573  K.  While 
it  is  heated  to  1573  K,  the  radius  decreases,  i.e.,  sample  warping 
increases.  This  level  of  warping  is,  however,  very  moderate  com¬ 
pared  with  that  observed  in  other  systems,  such  as  LTCC  [24]  and 
LSMZ/YSZ  [25].  This  is  mainly  due  to  the  fact  that  the  electrolyte 
GDC  layer  is  very  thin  compared  with  the  NiO-YSZ  anode  layer.  For 
a  thicker  GDC  layer,  the  large  difference  in  the  densification  kinetics 
between  the  two  materials  should  show  a  much  higher  curvature 
during  sintering. 


Fig.  11  plots  the  measured  rate  of  the  normalized  curvature 
change  as  a  function  of  the  isothermal  time  at  1573  K.  The  predic¬ 
tions  of  warpage  using  the  models  of  Cai  et  al.  [9]  and  Ranters  et  al. 
[15]  are  also  plotted  in  Fig.  11.  An  interesting  aspect  in  this  study 
is  that  the  same  experimental  data  were  used  for  both  models, 
allowing  a  direct  comparison  between  the  calculated  curves.  The 
measured  and  calculated  rates  are  in  good  agreement.  This  result 
indicates  that  continuum  mechanical  models  can  be  used  to  pre¬ 
dict  the  warpage  behaviour  of  a  solid  electrolyte  |  anode  bi-layer  if 
the  viscous  properties  of  the  components  are  properly  measured. 
The  model  of  Cai  et  al.  [9]  insignificantly  underestimates  the  camber 
rate  for  most  of  the  time  period,  whereas  the  camber  rate  calculated 
using  the  model  of  Ranters  et  al.  [  1 5  ]  is  larger  than  the  experimental 
value  for  the  first  200  s  of  isothermal  sintering.  An  apparently  better 
agreement  between  the  model  of  Cai  et  al.  [9]  and  the  measurement 
is  observed  in  this  study.  Previously,  Ravi  and  Green  [10]  measured 
viscoelastic  properties  by  cyclic  loading  and  calculated  the  warpage 
behaviour  of  laminated  alumina  tapes  with  different  green  densi¬ 
ties.  In  their  investigation,  however,  the  agreement  between  the 
experimental  and  theoretical  results  during  the  isothermal  sinter¬ 
ing  does  not  match  that  of  this  study.  On  the  other  hand,  Ranters 
et  al.  [15]  did  not  measure,  but  rather  estimated,  the  viscous  prop¬ 
erties  of  zirconia  layers  using  the  theoretical  model  of  Riedel  et  al. 
[26]. 
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Fig.  10.  Curvature  radius  of  bi-layer  measured  during  co-firing. 
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posite  is  determined  using  a  continuous  sinter  forge  technique  at  a 
temperature  of  1573  K. 

The  warpage  of  bi-layers  during  sintering  is  directly  observed 
in  a  furnace  equipped  with  a  CCD  camera.  During  co-firing,  the 
warpage  of  bi-layer  samples  with  a  thin  electrolyte  on  a  thick  bulk 
anode  is  moderate;  it  increases  at  the  beginning  and  decreases 
later.  The  warpage  behaviour  is  also  calculated  through  the  use 
of  the  measured  viscoelastic  properties  of  the  materials  for  two 
different  models.  The  calculated  results  are  in  good  quantitative 
agreement  with  the  observation,  in  particular  for  that  obtained 
using  the  model  of  Cai  et  al.  [9].  It  appears  that  the  accuracy  of 
constitutive  parameters,  in  particular  the  viscous  Poisson  ratio,  is 
critical  for  the  prediction  of  warpage  behaviour.  The  good  agree¬ 
ment  between  the  observed  and  predicted  warpage  in  the  present 
investigation  demonstrates  that  the  model  calculation  can  success¬ 
fully  be  applied  for  the  optimization  of  the  co-firing  of  bi-layer  SOFC 
components  when  the  viscous  properties  of  the  components  are 
properly  measured. 


A  number  of  reasons  can  be  considered  for  the  enhanced  agree¬ 
ment  in  the  present  study,  namely:  (i)  the  discontinuous  and 
continuous  sinter  forging  techniques  can  be  reliable  for  measuring 
the  sintering  parameters  for  high  density  and  intermediate  den¬ 
sity  regions,  respectively,  and  (ii)  the  camber  in  the  present  system 
is  much  lower  than  usually  observed  in  other  systems,  lowering 
the  discrepancy  between  the  model  and  the  experiment.  For  iden¬ 
tification  of  the  main  reason,  however,  further  experiments  with 
different  geometries  involving  other  systems  are  necessary. 

The  present  investigation  has  shown  that  the  warpage  behaviour 
of  an  asymmetric  half-cell  of  a  solid  electrolyte  and  an  anode  can 
be  predicted  well  through  the  use  of  adequate  continuum  mechan¬ 
ics  models  when  the  viscoelastic  properties  of  the  materials  are 
properly  measured  by  discontinuous  and  continuous  sinter  forging 
with  high  precision.  Once  the  sintering  and  viscoelastic  properties 
are  determined  by  a  sinter  forging  technique,  the  critical  design 
parameters  of  multi-layers  can  be  optimized  through  iterative  cal¬ 
culation  (temperature  program,  thickness  of  the  different  layers 
and  particle  size  of  initial  powder)  to  reduce  undesired  warpage. 
The  methodology  provided  in  the  current  work  is  also  expected  to 
be  useful  for  evaluations  and  control  of  the  warpage  behaviour  of 
other  multi-layer  structures  and  materials  used  in  SOFC  designs. 

5.  Conclusions 

The  viscoelastic  properties  of  gadolium-doped  ceria  (GDC) 
and  NiO-YSZ  are  determined  via  discontinuous  and  continuous 
sinter  forging  techniques  to  estimate  the  warpage  behaviour  of 
GDC|NiO-YSZ  bi-layers.  For  the  doped  ceria,  based  on  the  grain 
size  and  temperature  dependency  of  uniaxial  viscosity,  the  uniaxial 
viscosity  at  1573 1<  is  estimated.  The  viscous  Poisson  ratio  is  deter¬ 
mined  by  the  linear  relationship  between  the  applied  stress  and 
the  volumetric  densification.  The  ratio  estimated  using  the  model 
of  Venkatachari  and  Raj  [20]  is  found  to  be  in  good  agreement  with 
the  measured  value.  The  uniaxial  viscosity  of  the  NiO-YSZ  com- 
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